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ABSTRACT  s 

simulation  calculation  Is  carried  out  for  the  time  dependence  of  the  donor 
Intensity  as  one-step  acceptor  trapping  process  on  a  simple  cubic  lattice  of 
Interconnected  rods  of  dimensions  comparable  to  those  found  for  Itie  pores  of 
Vycor  (Corning  glass  7930).  The  donor  and  acceptor  are  allowed  to  occupy  random 
sites  on  the  surface  of  the  pore.  The  results-  are  fitted  using  the  equations  of 
Blumen  and  Klefter1  for  energy  transfer  on  a  fractal  structure.  An  "apparent 
fractal"  dimension,  <d’>,  is  determined  which  Is  shown  to  result  from  an 
excluded  volume  effect  and  Is  not  due  to  a  real  fractal  structure.  <d’>  Is 
found  to  depend  slightly  on  time,  the  acceptor  concentration,  and  the 
length  jwldth  ratio  of  the  pores. 

++  In  partial  fulfillment  of  the  Ph.  D.  degree 
4  on  leave  from:  E.R.A.  133<.N.R.S. , Uni  varsity  Paris  VI,  FRANCE 


INTRODUCTION 


Fractals^,  structures  with  a  dilation  symmetry,  have  attracted  a  great  deal 
of  attention  recently.  This  Is  due  to  their  utility  for  describing 

disordered  systems.  They  have  been  used  to  describe  systems  In 
polymer^*4,  protein4*®  and  surfaces  science®,  and  processes  such  as  crystal 
growth dielectric  breakdown®,  turbulence  and  chaos  .  Fractals  also  have  been 
used  to  describe  the  diffusion  of  liquids  Into  porous  media1®. 

11  12 

At  least  three  different  dimensions  are  required  to  define  a  fractal  *  . 

The  first  Is  the  Euclidean  dimension,  d.  In  which  the  structure  Is  embedded. 
The  second  is  called  the  fractal  dimension^,  D.  This  describes  the  dependence 
of  the  number  of  sites  N(R)  on  the  distance  R,  through  the  relation  (NCR)  *  R®). 
The  thlrc  dimension  is  the  spectral  or  fracton  dimension' ' ,  d,  which  governs 
the  random  walk  and  relaxation  processes  and  determines  the  density  of  states  of 
the  structure.  The  spectral  dimension  has  been  previously  discussed  in 
electron-spin  relaxation  studies  In  proteins4*®  and  triplet-triplet  annihilation 
studies  In  mixed  molecular  crystals1®*14.  More  recently,  studies  of  one  step 
electronic  energy  transfer  have  been  discussed,  both  theoretically1  and 

i  e 

exper Imentel  ly  ,  In  terms  of  the  fractal  dimension  D. 

Energy  transfer  between  Rhodamine  B  and  Malchite  green  doped  Into  Corning 

1 5 

Vycor  glass  (7930)  has  been  experimentally  analyzed  In  terms  of  the  Klefter 
and  Blumen1  equation  relating  the  fractal  dimension  to  the  temporal  behavior  of 
the  donor  Intensity.  A  fractal  dimension  of  I.74±0.12  has  been  determined. 
Electron  microscopic  studies1®  of  the  Corning  7930  Vycor  glass  showed  that  Vyoor 
glass  Is  a  rather  oompact  packing  of  spheres  whose  diameters  are  narrowly 
distributed  around  200  A.  This  packing  Is  consistent  with  other  results  such  as 


-the  pore  size  distribution,  which  is  found  to  be  sharp  (with  a  diameter  of 

t  «  c 

40  i  3  A)  ,  and  the  volume  fraction  of  pores,  which  constitutes  only  25-to-30  % 
of  the  total  volume1^. 

The  above  observations  together  with  the  fact  that  donor  and  acoeptor 

o 

molecules  have  finite  thicknesses  (/v5  A)  justify  representing  the  structure  of 

the  pores  In  Vycor  by  a  three  dimensional  simple  cubic  lattice  of  interconnected 

c  e. 

rods,  each  having  a  diameter  of  40  A  (the  pore  diameter)  and  a  length  of  200  A 

(the  diameter  of  the  spheres).  With  this  structure,  we  carried  out  simulation 

calculations  of  the  decay  of  the  donor  intensity  resulting  from  a  one  step 

dlpole-dl pole  energy  transfer  to  acceptors.  The  results  were  then  fitted  to  the 

Blumer.-Klefter  equation1  and  an  "apparent  fractal"  dimension  was  determined. 

The  latter,  as  expected  for  a  non-fractal  structure.  Is  found  to  depend  on  time, 

the  acceptor  concentration,  and  the  length:width  ratio.  However,  within  the 

range  of  times  and  concentrations  normally  available  to  experimentalist,  an 

"apparent  constant  fractal"  dimension  can  be  Inferred.  At  very  long  times 

and/or  at  very  lew  acceptor  concentrations  (conditions  that  cannot  be  reached 

experimentally  due  to  the  limitations  Imposed  by  the  value  of  the  radiative 

lifetime  of  the  donor),  a  dimensionality  of  3.0  Is  observed.  These  findings  are 

Ihe  result  of  the  excluded  volume  (of  the  spheres)  and  the  excluded  surface  (due 

to  surface  of  contact  between  spheres)  and  are  not  due  to  the  fractal  structure. 

The  Model 

We  consider  only  the  one  step  dipole-dipole  do  nor -to- accept  or  transfer. 
Donor-donor  energy  transfer  Is  excluded.  Furthermore,  the  energy  levels  of  the 
donor  end  acceptor  and  the  temperature  are  such  that  no  acoeptor  to  donor  back 
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transfer  I s  al lowed. 

Consider  the  case  of  one  excited  donor  located  at  a  position  Rp  which  Is 

—? 

surrounded  by  a  statistical  distribution  PA(Rp  ,  ra>  of  acceptors  at  position 

-*>  ^ 

Ra.  These  acceptors  can  trap  the  donor  excitation  In  a  trapping  time  6^(Ra-Rq) 

—+  — > 

which  depends  on  the  relative  distance  RA-RQ  between  the  acceptors  and  the 
donors.  It  can  be  shewn  that  the  decay  curve  of  the  donor  excitation  follow* the 
equation1  ^#1  ®: 


(t)  =  exp 


dlA 


,<V 


V 


[1  -  exp 


vVV 


(1) 


where  'aCQ  is  the  donor  lifetime  In  the  absence  of  acceptor  and  dRA  Is  the  volume 
element  over  which  the  Integration  Is  performed.  Equation  U)  has  been  used 
previously  assuming  all  donor  positions  are  equivalent.  In  our  structure,  Ihls 
Is  not  true  since  short  time  coupling  to  acceptor  distribution  Is  different  for 
donors  positioned  at  rod  Intersections  from,  those  situated  In  the  middle  of  the 
rods.  For  this  reason,  averaging  over  the  donor  position  deems  necessary.  The 
averaging  of  equation  (1)  over  the  positions  of  the  donors  with  distribution 


-> 


V° 


Fp(Rp)  leads  to  the  result: 

j  %  VV  “P  f  ♦  i  %  VV  RA>Ii-«xp  12) 

'  0  # 

.  !  R  .  -R__  1 

.  In  order  to  simplify  the  Integration,  we 


I  d*A  VV  RA)I1-CXP 

e  t  A  V 

'VS)1  r 

R 

o 

have  represented  the  narrow  distribution  of  the  radii  of  the  spheres,  by  a 


•  Ith  WV 


l 


regular  distribution  of  the  centers  of  the  spheres  on  a  simple  cubic  lattice, 
whose  spacing  Is  equal  to  twice  the  radius  of  a  sphere,  RQ.  Electronic  energy 
transfer  from  donors  to  acceptors  lying  either  on  the  same  sphere  or  on  two 

different  spheres  are  a  I  lowed.  As  no  molecule  can  occupy  sites  on  surface  of  a 

m 

contact  between  any  two  spheres  due  to  the  thickness  of  a  molecule  (**5  A),  an 
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excluded  surface  effect  Is  expected.  We  took  this  Into  account  by  transforming 
the  lattice  of  spheres  Into  a  lattice  of  rods.  Interconnected  and  directed  along 
the  three  axis  of  the  cubic  lattice.  The  lengths  of  the  rods  are  equal  to  the 
diameter  of  the  sphere  and  their  d! antlers  are  equal  to  the  pore  diameters. 
Donors  and  acceptors  are  allowed  to  occupy  random  sites  on  the  surface  of  these 
rods. 


RESULTS  AND  DISCUSSION 


The  simulation  consists  of  calculating  lp(t)  e(+^f  *  using  Eq.(2)  at 
different  times  t;  for  different  acceptor  concentrations,  p^;  and  for 


different  rod  length :wldth  ratio  (L/W).  The  numerical  Integration  was  carried 
out  on  the  surface  where  the  donors  and  acceptors  are  randomly  distributed. 
This  represents  the  time  dependence  of  the  donor  Intensity  resulting  only  from! 
energy  transfer  processes  (and  not  from  first-order  decay). 

In  order  to  fit  our  simulation,  we  have  used  the  equation  previously 


derived  : 


1  (t)  -  exp  JL  +  p  a  (— ) 
V  .A  1  _ 


t  .d  (D) / s 


(») 


.i 


For  dipole-dipole  transfer,  s  »  6;  A  Is  a  time  independent  function;  and  d  Is 
the  Euclidean  or  fractal  (D)  dimensionality  depending  on  the  structure  of  the 
system  studied.  This  equation  has  been  used  by  Even  et  el1^  to  fit  experimental 
data  for  energy  transfer  between  Rhodamlne  B  and  Malctiite  green  on  Vycor  and 
gave  D  •  1 .74  ♦  0.12. 

It  is  obvious  that  the  simulation  of  the  temporal  behavior  of  the  donor 
intensity  tr.  our  model  should  not,  and  Indeed  does  not,  follow  a  simple  decay 
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law  of  the  form  of  Eq.  (3)  In  a  strict  sense.  Thus,  the  exponent  of  the  time 
In  Eq.  (3)  cannot  be  D,  since  we  selected  a  cubic  and  not  a  fractal  lattice. 
As  w  1 1 1  be  shown  below,  this  exponent  does  not  have  to  be  Integer  at  all  times 
or  acceptor  concentrations;  thus.  It  cannot  be  d.  Let  us  call  It  d' .  It  Is 
found,  as  expected,  that  d*  depends  on  the  time  and  on  the  acceptor 

concentration.  In  order  to  determine  the  time  and  concentration  dependences  of 
th  is  exponent,  the  negative  logarithm  of  the  calculated  Intensity  has  been 
plotted  versus  time  on  a  In- In  plot  for  different  acceptor  concentrations  (see 
fig.  1).  In  such  plots,  the  slopes  of  the  curves  are  equal  to  d'/6.  Although 
Fig.  1  shows  that  d'/6  depends  on  the  time.  It  also  shows  that  It  Is 

approximately  constant  for  a  given  range  of  time,  concentration,  and  L:W.  This 
allows  us  to  use  Eq.  (3)  to  determine  an  average  value  of  d' ,  <d'>,  for 
structures  such  as  the  one  we  are  presently  using. 

It  is  then  possible  that,  on  Vycor-llke  structures,  the  decay  can  be 

described  with  a  pseudo-constant  exponent.  For  Instance,  nost  experimental 
results  are  obtained  during  the  time  required  for  the  Initial  Intensity  to 

decrease  by  at  most  a  factor  of  100.  During  this  time  range,  the  value  <d’>  is 
apparently  constant. 

From  Fig.  (1),  the  best  value  of  <d’>  Is  determined  for  the  different  time 
and  concentration  range  of  interest.  With  these  values,  equation  (3)  Is  then 
tested;  and  the  results  are  shown  In  Fig.  (2)  which  gives  a  straight  line  for 
In  I  vs.  t<<3’>/6. 

Table  1  gives  the  values  of  <d’>  for  the  different  concentrations  and  LtW 
for  the  time  range  required  for  the  donor  Intensity  to  decay  to  1/100  of  Its 
Initial  value.  In  this  table,  the  transfer  decay  times  In  units  of  a 

critical  transfer  time  (^>,  the  time  of  transfer  at  a  do  nor -accept  or  distance 


of  the  sphere  radius,  Rc,  are  also  given  for  different  concentrations  and  LsW 
values.  As  expected,  as  the  L:W  ratio  Increases,  <d*>  decreases,  and, 
eventually,  one-d  Intensions  I  behavior  is  observed.  Also,  as  the  concentration 
decreases  (or  the  tine  Increases),  the  effects  of  the  excluded  volume  and 
surface  disappears;  and  three-dimensional  behavior  Is  recovered. 

The  approximate  simulation  above  shows  that  <d'>  *  1.7  (the  value 
determined  experimentally)1^  can  be  obtained  for  an  acceptor  concentration  of  6 
molecules  per  sphere  and  an  L:K  ratio  of  10.  These  values  have  to  be  compared 
to  the  experimental  ones:  If  all  of  Ihe  pores  between  the  spheres  In  Vycor  are 
Interconnected,  the  L;W  ratio  should  be  200:40  «  5.  The  concentration  used 
by  Even  et  al1^  was  2-to-3  molecules/sphere.  If  one  assumes  a  random 
distribution  of  donors  and  acceptors  within  the  Vyoor  glass. 

Our  model  oversimplifies  the  Vycor  structure.  The  silicate  spheres 
probably  are  randomly  packed,  and  the  microstructures  of  the  pores  ere  not 
Known.  Thus,  one  should  not  take  our  comparison  with  experiments  too  seriously. 
It  Is  not  the  purpose  of  this  study,  however,  to  attempt  to  describe  the  Vycor 
structure,  nor  to  compare  our  results  with  experimental  observation.  The  main 
purpose  of  this  work  Is  to  show  that.  If  one  studies  one-step  trapping  process 
on  a  non-fractal  (e.g.,  cubic)  lattice  with  excluded  volumes  of  length  scales 
equal  or  greater  than  the  characteristic  energy  transfer  length,  an  apparent 
fractal  behavior  can  be  obtained  over  limited  range  of  time  (or  acceptor 
concentrations).  Of  course,  one  can  try  tc  distinguish  between  a  true  fractal 
and  an  apparent  fractal  behavior  If  the  experimental  study  is  carried  out  over  a 
long  time  scale.  Most  unfortunately,  this  Is  not  possible  due  to  the  limitation 
Imposed  on  the  experimentalist  by  the  finite  lifetime  of  the  excited  state  of 


the  donor 


Extensions  of  Ihls  work  to  simulations 


of  other  Interesting 


structures  are  no*  In  progress,  and  the  details  of  these  si  mu!  attorn*  1 1 1  be 
published  later. 
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Figure  caption: 

Fig.  1: 

Plot  of  the  In(-lndntenslty))  versus  In(tlme)  for  different  acceptor 
concentration  on  a  simple  cubic  lattice  of  rods  of  length:* Idth  ratio  of  5 
for  top  and  10  for  the  bottom  figure.  Dipolar  interaction  (s*6)  Is  assumed.  Ip 
denotes  the  donor  luminescence  Intensity;  the  slope  of  the  curves  defines  an 
exponent  d'/s  which  varies  with  time,  and  concentration  as  expected  for 

non- fractal  structure.  PA  •  20,  PA  «  6,  PA  -  2.  t  Is  In 

arbitrary  unit. 

Fig.  2: 

Plots  of  the  calculated  In  of  the  donor  Intensity  versus  t<<J,>^6  In  a  semi-log 
scale,  there  <d'>  Is  the  average  fractional  dimension  value  obtained  from  the 
averaging  slope  of  Fig.  1.  These  results  show  that  "fractal- 1  Ike n  behaviors  can 
be  observed  over  the  range  of  time  and  concentration  used  experimentally  If  a 
structure  similar  to  the  one  with  which  we  approximated  the  Vycor  structure  Is 
used.  *  20,  P.  ■  6,  P.  ■  2.  t  Is  In  arbitrary  unit. 
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Table  1:  The  dependence  of  <d*>  (the  average  pseudo  dimensionality)  and  the 


transfer  time  on  the  acceptor  concentration  PA  and  length:wldth  ratio  L/W. 


Is  the  time  for  which  the  donor  Intensity  has  decayed  to  1/e  of  Its  Initial 


value  In  unit  of  ^ ,  the  characteristic  time  for  a  transfer  between  a  donor 


and  an  acceptor  separated  by  R  the  radius  of  the  sphere. 


<d'» 


2.1  2 


.6  10 


2.15  1.7  1.6 
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